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A lattice model for a symmetrical copolymer "bottle brush" molecule, where two types (A,B) of 
flexible side chains are grafted with one chain end to a rigid backbone, is studied by a variant of the 
pruned-enriched Rosenbluth method (PERM), allowing for simultaneous growth of all side chains 
in the Monte Carlo sampling. Choosing repulsive binary interactions between unlike monomers and 
varying the solvent quality, it is found that phase separation into an A-rich part of the cylindrical 
molecule and a B-rich part can occur only locally. Long range order (in the direction of the backbone) 
does not occur, and hence the transition from the randomly mixed state of the bottle brush to the 
phase separated structure is strongly rounded, in contrast to corresponding mean field predictions. 
This lack of a phase transition can be understood from an analogy with spin models in one space 
dimension. We predict that the range of microphase separation along the bottle brush backbone 
can be controlled on the nanoscale by varying the solvent quality. 



I. INTRODUCTION AND MOTIVATION 

Macromolecules with a comb-like architecture, where 
side chains are grafted to a linear (backbone) chain, find 
increasing interest: if the backbone is a rigid polymer, its 
solubility and processibility is improved [J M y] ; if the 
grafting density of the side chains is very high, cylindrical 
"bottle brush" -shaped objects form, which under certain 
conditions show a thermally induced collapse transition 
to a spherical structure, providing interesting perspec- 
tives for the design of " molecular actuators" [3] ; if copoly- 
mer bottle brushes are created by dense grafting of two 
types of side chains (A, B) on the backbone, interesting 
"horseshoe" and " meander" -like structures are observed, 
which is attributed to a local phase separation along the 
backbone of the bottle brush Clearly, this creates 
additional possibilities when one considers the desired [|[ 
use of such shape persistent cylindrical macromolecules 
as building blocks in molecular assemblies. Thus, the 
problem of phase separation within a simple copolymer 
bottle brush has also found attention both by analytical 
theory [f| and computer simulation Q > the possibility 
to create "Janus cylinders" (upper half of the cylinder 
containing the B monomers, lower half containing the A 
monomers) was suggested, and predictions how the phase 
transition point from the randomly mixed state (where 
both A and B monomers are homogeneously distributed 
in the cylinder volume) to the separated state depends 
on the chain length of the side chains N were given @. 

In the case of flexible backbones there is a delicate in- 
terplay between the local segregation between the A and 
B side chains and the global geometric structure (such 
as "horseshoe", etc.) of the bottle brush, requiring a 
detailed and accurate understanding of the phase sepa- 
ration in these molecules. We reconsider this problem, 
performing_Monte Carlo simulations of a similar model 
as in Ref. [7(, but our analysis focuses on a different as- 
pect, that has been completely ignored in the literature 
so far: single bottle brush molecules with side chains 
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FIG. I: Crossover scaling behavior of the correlation length £ 
as a function of temperature T. £ cx {T/T c - l)~ u for T > T c , 
and £ cx R 2 F(T) /k B T for T < T c . 



of finite chain length N and very long, ideally infinitely 
long length of the backbone chain are quasi-one dimen- 
sional (Id) objects. Statistical mechanics of systems with 
short range interactions implies quite generally [1, H[ that 
there cannot be any sharp phase transitions in Id sys- 
tems in thermal equilibrium. A point in case is e.g. the 
Id XY model, a chain of spins where each spin i is de- 
scribed by an angle (fi in the xy-plane, with < ifi < 2ir. 
Neighboring spins along the chain experience a coupling 
— Jcos((fi — (fj). While mean field theories predict that 
ferromagnetic order along the chain occurs for temper- 
atures T < T C MF = cJ, where c is a constant of order 
unity [H , an exact solution of this problem [H, Q shows 
that ferromagnetic long range order is unstable against 
long-wavelength fluctuations, and actually the ferromag- 
netic correlation length £ grows completely gradually as 
the temperature is lowered, 

£ = 2a{J/k B T) , (1) 
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a being the lattice spacing. Thus T c — 0: note that £ 
has to diverge when T c is approached, and the singular- 
ity of £ predicted by mean field theory at T = T^ tF 
is completely washed out. This consideration can be 
generalized to cylinders of cross section area nR 2 [Io| . 
Then Eq. |T]) gets replaced by a more complicated be- 
havior that is sketched in Fig. [T] For T > T c (the critical 
temperature of the corresponding bulk three-dimensional 
system) £ oc (T/T c — 1) ^, where v is the appropriate 
critical exponent, until £ is of the order of R. Then the 
critical singularity is rounded off, and a crossover sets in 
to a relation analogous to Eq. |T]) 

£ = 2T(T)irR 2 /k B T oc {R 2 /a){J/k B T) , T^O. (2) 

Note that the "spinwave stiffness" (helicity modules) 
r(T) characterizes the cost of long wavelength order pa- 
rameter rotations, and T(T — > 0) oc J. In the present 
work we shall provide evidence that the correlation length 
describing phase separation of a "Janus cylinder" type 
behaves qualitatively similar to Eqs. {TJ and @, cf. 
fig.[2p. Then the Flory-Huggins parameter xab that de- 
scribes the incompatibility between A and B monomers 
takes the role that (J/k B T) plays for the spin model. 
However, the problem whether for N — > oo (and depend- 
ing on the solvent condition) sharp phase transitions are 
restored is nontrivial, and in particular, describing the 
extent of rounding of the (mean field) transitions de- 
scribed in [6(] by these long range fluctuations along the 
chain backbone remains a challenge. 

II. MODEL AND SIMULATION TECHNIQUE 

As considered previously [f| 0] , the generic model for 
this problem treats a completely rigid rod, oriented along 
the z-axis of a simple cubic lattice, and the side chains 
of type A and B are grafted in an alternating way, every 
site of the backbone chain being the grafting site of a 
side chain end. Of course, creating a statistical copoly- 
mer where the type of side chain along the backbone is 
chosen at random would be more realistic, from the ex- 
perimental point of view. However, the quenched disor- 
der thus introduced would encumber a theoretical treat- 
ment considerably, therefore this has not been considered 
previously @, 0] , and hence we study here regularly alter- 
nating copolymers, for the sake of better comparison to 
previous work. Also, we introduce a periodic boundary 
condition in z-direction, to avoid end effects associated 
with a finite backbone length. 

The side chains are self-and mutually avoiding ran- 
dom walks. Between monomers of the side chains and 
monomers of the backbone chain only the excluded vol- 
ume effect is considered. We allow for nearest neighbor 
interactions e AB , £aa — £bb = e between the respec- 
tive pairs of monomers of the side chains, and hence the 
partition sum for this model is 

Z = J2<l mAA+mBB <lAB B > ( 3 ) 




FIG. 2: Perfect phase separation of side chains in a binary 
(A, B) copolymer bottle brush with alternating grafting se- 
quence ABAB ■ ■ ■ into a " Janus cylinder" structure implies 
formation of an interface plane (part a), shaded) between the 
^4-rich part (bottom) and the B-rich part (top) of the cylin- 
drical bottle brush. The local orientation of this interface can 
be characterized by a vector oriented normal to it (arrows). 
At nonzero but low temperature phase separation will occur 
locally, but entropy will lead to long-wavelength fluctuations 
of the orientation of this vector (b), destroying axial long 
range order along the direction of the backbone of the bottle 
brush. Since the interface is not a sharp plane at nonzero 
temperature, but rather has a finite width, a numerical char- 
acterization of the local orientation of the interface is difficult. 
An essentially equivalent but numerically unambiguous char- 
acterization of this "Janus cylinder" -type ordering is achieved 
by calculating the unit vector from the grafting site (n) of each 
chain to the projection of its center of mass position into the 
xy- plane R" m (n), see part (c). 



where q = exp(-e/fc B T), q AB = cxp(-e AB /k B T), and 
m AA , mBB, m AB are the numbers of non-bonded occu- 
pied nearest neighbor monomer pairs A A, BB, and AB, 
respectively. The sum in Eq. §5§ extends over all possi- 
ble configurations of these walks. Varying q in the range 
1 < q < 1.5, we cover the full range from good solvents 
(e = 0) to poor solvents (e = 1.5), since the ^-solvent cor- 
responds to a choice [ll( qg = exp(— e/k B To) ~ 1.3087. 
We vary q AB in the range < q A B < q- The case q A B = 
corresponds to a very strong repulsion between A and B, 
while for q AB = q the chemical incompatibility vanishes 
{recall that xab oc e AB - (e AA + e BB )/2}. 

For our simulations, we use the pruned-enriched 
Rosenbluth model (PERM) [ill- This is a biased 
chain growth algorithm with population control. PERM 
has been applied very successfully both to linear and 
branched polymers (stars [H, lattice animals [HI]). 
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FIG. 3: Correlation function C n describing the correlation of local phase separation along the backbone plotted vs. the distance 
n on a semi-log scale, for good solvent conditions (e = 0), N = 6 (left part) and N = 18 (right part). All data are for Lb = 32, 
and several choices of qAB are included, as indicated. Note that C n = (< ■ Sf +n > + < SP +1 ■ SP +1+n >)/2. Curves show 
the best fit to the data, i.e. C n ~ [exp(— n/£) +exp( — (L& — n)/£)] 




Similar to [H, [HI, the bottle brush is generated by 
adding one monomer to each side chain until all side 
chains have the same number of monomers, thus growing 
all side chains simultaneously. For efficiency, side chains 
are grown with higher probabilities in the directions per- 
pendicular to the backbone and in the direction where 
there are more free next neighbor sites. This additional 
bias must be taken into account by suitable weight fac- 
tors. High statistics (10 5 ~ 10 6 independent configura- 
tions) is obtained by using the population control (prun- 
ing/cloning) [ll|, [lj, 0, EH- The backbone length also 
is varied (L& = 32, 48, 64) to check that finite size effects 
caused by too small are still negligible. 



III. RESULTS AND DISCUSSIONS 

As has also been discussed in [7[ , the quantification of 
the phase separation that occurs in a bottle brush is diffi- 
cult. We follow here an approach that is motivated by the 



analogy with the spin problem mentioned in the introduc- 
tion (fig. ^ . We draw from each grafting site the vector 
R" m to the center of mass (CM) of each side chain of 
type a (a = A or B), and use the x and y-components of 
these vectors to define a unit vector S!j* that points along 
that direction in the xy-plane. Making use of the fact 
that we have grafted A and B chains in an alternating 
way, we compute the correlation functions < Sf • S^_ n > , 

< Sf +1 ■ S? +1+n >, for n = 0, 2, 4, 6, .... If we would 
have a symmetry breaking of "Janus cylinder" type (i.e. 
a plane containing the backbone exists, such that the 
A-monomers are on one side of this plane and the B- 
monomers on the other), the above correlations would 
not decay to zero as n increases, but rather settle down 
at a nonzero order parameter square. Figs.[3]and Ushow 
that this is not the case: for good solvent conditions q = 1 
even strong incompatibility (qAB ~^ 0) can induce only a 
correlation of rather short range. In poor solvents there 
is also little correlation when the side chains arc very 
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FIG. 5: Snapshots of the copolymer bottle brushes for a poor solvent (q = 1.5), L;, = 64, N = 18, qAB = 1.5 (mixed state) 
and qAB = 0.1 (state with local phase separation). Monomers A , monomers B, and monomers on the backbone are shown in 
black, gray, and white colors, respectively. 
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FIG. 6: Inverse of correlation lengths £ vs. Xab f° r copolymer bottle brushes with backbone length Lb = 32, 48, and 64, and 
three chain lengths N = 6, 12, and 18. Both good (left) and poor (right) solvent conditions are shown. 



short (N — 6) and much stronger correlations occur for 
longer chains (N = 18). Although our data are already 
somewhat affected by the finite length of the backbone 
(note the symmetry of C n = CL h - n ), they are fitted by 
exp(— n/£) + cxp(— (Lb — ?i)/£) quite well. In fig. corre- 
sponding snapshots (similar to those of 0]) show rather 
pronounced phase separation, but this clearly is an arti- 
fact of the finite system size. The correlation lengths £ 
which are determined by fitting C n for n > 2 as shown 
in fig. [3] and [4] are shown in fig. [6l As expected in view 
of the discussion given in the introduction, £ increases 
gradually with increasing chemical incompatibility xab 
(decreasing temperature). Choosing poor solvent and/or 
N larger makes £ also larger. Since l/£ decreases linearly 
as Xab ~ * 0, one can extrapolate the data to Xab = 
by fitting a straight line l/£ = 1/£at + b/xAB- Results 
of £n in a poor solvent are shown in fig. [7] It gives an 
indication that a sharp phase transition (l/£/v - ► 0) will 
develop in the limit of large xab and large N. For very 
poor solvent condition, one expects that the brush fills 
a cylinder of radius proportional to y/~N densely, so the 
situation should be analogous to Eq. ((2]), with irR 2 being 
replaced by N. Our data are roughly compatible with 
such an interpretation. However, the radial density pro- 
files of the bottle brush molecules (fig. [7|) show that the 
actual density profile p(r) is still far from a profile of a 
densely filled cylinder (p(r) = 1 for < r < R). This low 
density is probably responsible for the fact that £ stays 
finite as Xab ~ * (fig- EJ) . 

Alternatively, one could partition a bottle brush along 



its backbone direction into disks with height a (distance 
between the neighboring grafting sites in our model). 
Each disk only contains one grafting site. The unit vec- 
tor Sf is then defined by the direction from each grafting 
site to the CM of monomers of type a (a = A or B) in 
the x-y plane. However, our results show that the corre- 
lation functions of this observable decay much faster due 
to the large fluctuations of the orientation of the vector 
Sf. 

Summarizing these results, we have presented evidence 
that the phase separation towards a "Janus cylinder" 
structure develops completely gradually as the incom- 
patibility between the two types of monomers increases. 
For our short side chain lengths, no trace of the sharp 
phase transitions predicted by Stepanyan et al. @ could 
be detected (a similar conclusion was also reached in 0]). 
Very large length of side chains and/or very strong incom- 
patibility is needed to produce a correlation extending 
over a large number of side chains. The fluctuation that 
has been ignored by the mean field theory [|| , long wave 
length random "twist" -like rotation of the local inter- 
face between the j4-rich and B-rich region in the " Janus 
cylinder" , costs very little energy and destroys long range 
order dramatically (fig. Wp)- On the other hand, varying 
suitably side chain length and solvent quality one can 
control the nanoscopic length scale over which Janus-type 
phase separation occurs along the backbone of the bottle 
brush. For asymmetric side chain lengths (and flexible 
backbones) this length scale is related to the spontaneous 
curvature of the cylinder-shaped bottle brush polymer. 
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FIG. 7: Inverse of correlation length £jv in the limit of xab — » oo vs. 1/N for a poor solvent (<j = 1.5). The straight line 
shows the asymptotic behavior as N — > oo (left part). Monomer density profile p(r) plotted vs. r for a poor solvent (g = 1.5), 
Qab = 0.1 and Lb = 64 (right part). 



Controlling this curvature is desirable for the anticipated 
application of such molecules in nanotechnology [5j . 
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